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Abstract. We elaborate on a recently conjectured relation of Painleve transcendents 
and 2D CFT. General solutions of Painleve VI, V and III are expressed in terms of 
c = 1 conformal blocks and their irregular limits, AGT-related to instanton partition 
functions in J\f = 2 supersymmetric gauge theories with Nf — 0,1,2,3,4. Resulting 
combinatorial series representations of Painleve functions provide an efficient tool for 
their numerical computation at finite values of the argument. The series involve 
sums over bipartitions which in the simplest cases coincide with Gessel expansions of 
certain Toeplitz determinants. Considered applications include Fredholm determinants 
of classical integrable kernels, scaled gap probability in the bulk of the GUE, and all-order 
conformal perturbation theory expansions of correlation functions in the sine-Gordon 
field theory at the free-fermion point. 



1. Introduction 

Painleve transcendents [13] are nowadays widely recognized as important special functions 
with a broad range of applications including integrable models, combinatorics and 
random matrix theory. Many aspects of Painleve equations, such as their analytic and 
geometric properties, asymptotic problems, special solutions and discretization, have been 
extensively studied in the last four decades. 

From the point of view of the theory of classical special functions (58] , the surprising 
feature of these developments is the absence of transparent connection to representation 
theory. Instead, the Riemann-Hilbert approach [22] is typically used. It is well-known 
that Painleve equations emerge most naturally in the study of monodromy preserving 
deformations of linear ODEs. Thus, by analogy with the solution of classical integrable 
systems by the inverse scattering method, the questions on nonlinear Painleve functions 
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may be asked in terms of linear monodromy. In particular, one may attempt to realize 
the following program: 

• label different Painleve functions by monodromy data of the auxiliary linear problem, 

• express their asymptotics near the critical points in terms of monodromy, 

• construct full solution using the asymptotic behaviour as initial condition. 

Starting from the foundational work of Jimbo |31| . there are many results available 
on the first two points, but the lack of algebraic structure makes the last one difficult to 
tackle. In other words, the question is 

...how does one combine asymptotic information about the solutions obtained from the 
Riemann- Hilbert problem, together with efficient numerical codes in order to compute 
the solution u(x) at finite values of x? [14, Painleve Project Problem]. 

In [27], a solution of this problem was suggested for the sixth Painleve equation. 
It was shown that Painleve VI tau function r VI (i) can be thought of as a correlation 
function of primary fields in 2D conformal field theory [5] with central charge c = 1. 
Under natural minimal assumptions on primary content of the theory and fusion rules, 
t vi (t) may then be written as a linear combination of Virasoro conformal blocks. Being 
purely representation-theoretic quantities, these CFT special functions can be computed 
in several ways. In particular, the recently proven pQ AGT conjecture [2] relates them 
to instanton partition functions in M = 2 SUSY 4D Yang- Mills theories (6j EU H3j E] , 
which can be expressed as sums over tuples of partitions. This results into combinatorial 
series representations of r VI (t) around the critical points 0, 1, oo. 

The aim of this note is to extend the results of (27] to Painleve V and Painleve III, and 
to make them accessible to a wider audience. With this purpose in mind, we deliberately 
include some background material and illustrate our claims with a number of explicit 
examples and applications to random matrix theory and integrable QFT. 

The plan is as follows. Section [2] sets the notation and explains the relation between 
different Painleve equations and their various forms. In Section |3j we recall some 
basics on conformal blocks and AGT correspondence. Conjectural general solutions of 
Painleve VI, V and III are presented and discussed in Section |4| In particular, it is 
shown that our combinatorial expansions can be seen as a generalization of Gessel's 
theorem representation of classical Toeplitz determinant solutions. Section [5] is devoted 
to applications, which include Fredholm determinants of classical integrable kernels 
(hypergeometric, Whittaker, confluent hypergeometric and modified Bessel), scaled GUE 
bulk gap probability and correlators of exponential fields in the sine-Gordon model at the 
free-fermion point. 
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2. Painleve equations 

2. 1 . Conventional form 

Painleve VI, V, and III (P VI , P v , P m ) first appeared as a part of the classification of 2nd 
order, 1st degree nonlinear ODEs without movable critical points. In this context, they 
are usually written as follows: 

Painleve VI: 



dt 2 2 \q q-1 q - t J \dt J \t t-1 q - t J dt 
, 2q(q-l)(q-t) ( /3t 7 (t - 1) St(t - 1) 



t 2 (t — l) 2 V q 2 (q-1) 2 (q-t) 2 
Painleve V: 

d 2 q (1 1 \(dqV ldq , [q-1) 2 ( 8\ 1Q 5q(q + 1) 

d^ = {Yq + q—l){-d-t) --tTt + —^ 2 —\ aq+ -q) + T + ^l-' ^ 
Painleve III: 

d 2 q 1 ( dq\ 2 ldq aq 2 + 8 , 5 , 

1 ] 1 + * +7? + -• (2.3) 



dt 2 q \dt J t dt t q 

It is often convenient to use instead of P m an equivalent equation, 
Painleve IIP: 

d 2 g = 1 fdqV ldq q 2 (a + ^q) 8 5 

dt 2 q\dt) tdt At 2 At Aq { ' } 

which reduces to P m by setting t m , = t 2 n , q m , = t m q m . 



2.2. Parameterization 

We write four P VI parameters as 

(«, 0, 7, S) yi = ( ((900 + i) 2 , -# 2 , # 2 , - £ 2 ) . (2.5) 

If 5 7^ in P v , then one can set 5 = —\ by rescaling the independent variable. P v 
with 5 = is reducible to P UI (see e.g. transformations (1.24)-(1.26) in [12]) which will 
be treated separately. Hence we may set 

(a,/3, 7 ,5) v = (20jj,-20 t a ,20*-l,-i). (2.6) 

The case of P m is slightly more involved. In the generic situation, when 7$ 7^ 0, one 
can assume that 7 = — 5 = A by rescaling t and g. We will then write 

(a, /3,7,5V =(8^,4-8^,4,-4). (2.7) 
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The variable change q — > q 1 maps P m with 5 = to P m with 7 = 0. Assume that 7 = 
and a>5 ^ 0, then the scaling freedom can be used to set 

(a, /?, 7, 8) m2 = (8, 4-80,, 0,-4). (2.8) 

For 7 = 5 = 0, a(3 ^ we can set 

(a,/3, 7 ,5) IIl3 = (8,-8,0,0). (2.9) 

Finally, for a = 7 = (and, similarly, for (3 = 5 = by q— > q~ x ), the general (two- 
parameter) solution of P m is known [SB]- It reads 

q(t)=fit +^t + — — t + , 



where /1, z/ are two arbitrary integration constants. Excluding this last solvable case, there 
remain three inequivalent -P m 's with two, one and zero parameters. Significance of the 



degenerate equations P ul2 and P m was realized in (49] from a geometric viewpoint, and 
later they were extensively studied in |46j. 

2. 3. Hamiltonian form 

Painleve equations can be written as non-autonomous hamiltonian systems |37J. In 



this approach, (2.1), (2.2) and (2.4) are obtained by eliminating momentum p from the 
equations 

dq_dH 1 dp__dll ± J -VI VIII' 

dt ~ dp ' dt~ dg> J- vi,v,m li2)3 , 

where the Hamiltonians are given by 

20 o 20i 20 t -l 



t(t-l)H vl =q(q-l)(q-t)p(p- 

\ Q q~ 1 q-t 

+ (0 O + t + 0i + 0oo) (0 O + 6 t + 0i - 0oo - 1) g, (2.10) 
= (g - l)(pg - 20 t )(pq -p + 20,) - tpq + ((0, + t ) 2 - 2 ) q + 

+ (0 t -|)t-2(0, + |) 2 , (2.11) 

tH u , = {pq + 0,) 2 + tp - 0^ - q -, (2.12) 
t# m , = (pq + 9^ 2 + tp-q, (2.13) 

^ m ,=pV-?-^. (2-14) 

The hamiltonian structure is crucial for the construction of Okamoto-Backlund transfor- 
mations [451 . generating an infinite number of Painleve solutions from a given one. 
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The time-dependent Hamiltonians (2.10)-(2.14) themselves satisfy nonlinear 2nd order 
ODEs. To write them, introduce auxiliary functions 

a VI = t(t - 1)H V1 - q(q - l)p + (9 + 6 t + 9 1 + 9^) q 



fa , g ^2 , , #1 + C - - A? ~ 
-(^O + ^lJ t-\ , 

a 3 = tH 3 , J = V,IIIi A3 . 

They satisfy the so-called a- form of Painleve equations [231 ED] : 



P, 



a'[t{t-l)a") + 2a'(ta'-a) 



a 



/\2 



1 2 



(2.15) 
(2.16) 

(2.17) 



P 



(a' + (9 t + Oj 2 ) (a' + (9 t - e^f) (a' + (9 + 9 x f) (a' + (9 - 9 1 ) 2 ) 

2 1 



P 



(ta") 2 = (a-ta' + 2 (a') 2 ) - - {(2o-' - 9*f- A9 2 ^ {(2a 1 + 0*) 2 - A9 2 } , (2.18) 

(ta") 2 = (4 (a') 2 - 1) (a - ta') - 40Aa' + (9 2 + 0*) , ( 2 .19) 

(ta") 2 = 4 (a'f (a - ta') - 49*a' + 1, (2.20) 

'"2 

P m , : (ta") 2 = 4 (a') 2 (a - ta') - 4a', (2.21) 

which also appear in the classification of 2nd order, 2nd degree ODEs with Painleve 
property 



The solutions of (2.1)-(2.4) can thus be mapped to solutions of (2.17)-(2.21 ). 
Conversely, one can recover conventional Painleve functions from the solutions of a- 
Painleve equations using the following formulas: 



P- 



1 / 1 1 

+ 



q-t 2 \t t-l 
_ 2^t(t - l)q" + (a' + 9 2 - gfj ((2t - 1) a' - 2a + 9j - g) + 41%, (gg - 9\ 
2t(t - 1) (a' + (9 t - O^f) {a' + (9 t + ^) 2 ) 
2(ta" + a-ta' + 2(a') 2 ) 



(2.22) 



P 



P 



q 

Q 
Q 
Q 



(2a' -9*)- A9 2 
2ta" + 49*a' - 29* 

4(a') 2 - 1 ' 
ta" + 29, a' - 1 



2 (a 



a' 



Finally, define the tau functions of P VI , P v and P by 



a VI (t)= t(t -1)^ In [t 



/i2 I /j2 /i2 /i2 



/}2_i_/}2 /i2 /)2 

(1-*) 2 ^(i) 



(2.23) 
(2.24) 
(2.25) 
(2.26) 

(2.27) 
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,(t)= t|ln(e-^t-^ 2 -*. v (t)) 



(2.28) 



J = III' 1)2i3 . (2.29) 

Our solution below is formulated in terms of combinatorial expansions of these tau 
functions in powers of t. Expansions of cr's and q's can then be obtained from the relations 



(2.22)-(2.26) and (2.27)-(2.29) 



2. 5. Coalescence 

As is well-known, Painleve VI produces all other Painleve equations in certain scaling 
limits. The equations considered in the present paper form the first line of the coalescence 
cascade 

VI > V > IIIi > III 2 — > III3 



IV 



II 



I 



Fig. 1: Coalescence diagram for Painleve equations 

Every step to the right or to the bottom of the diagram decreases by 1 the number of 
parameters in the corresponding equation. 

Let us now describe the scaling limits we need (1st line) and the transition P m — > P x 
which seems to be missing in the literature (cf. e.g. the degeneration scheme in |46j). 



p vi -> p v - set in P VI 



0i 



A + 6, 



A -0* 



2 ' 2 ' 

then solutions of P„ can be obtained as the limits 



1 - q v (t) 

°v (t) 



t/A 



lim 

A— >oo 



lim 

Qyn (*/ A ) 
A 2 



(j VI (t/A) 



r v (t) = lim (t/Af +e * 

A— >oo 



T, . 



4 

(*/A) , 



V in' 



this limiting transition is described by 
A - 6+ . A + 0* 



0o 

^III' (t) 



<z v (t/A)) 



0t 

lim A (1 - 

A— >oo 

/ A 2 + ^ + 2 

hm _^ 1 + av (t/A) 

A— >oo \ z 



V (*) = A lim T v (*/A) 



(2.30) 

(2.31) 

(2.32) 
(2.33) 

(2.34) 
(2.35) 

(2.36) 
(2.37) 
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q **>® = & e *V (*/**), (2 - 38) 

^(*) = ^ V ( * /6, * ) ' (2 ' 39) 

r ni/ (t) = lim r m/ (t/0,) . (2.40) 

m P . 

III' 2 Ill'g 

^ = , h i^ (t/ ^ ) ' (2 ' 4l) 

r m , (t) = lim r m , (t/0,) . (2.43) 

0* = 3At, ^ 2 =16A-(l + ^), (2.44) 

^ (^) = 2A ^) + 8A§ + \ A ~Kv ( 2 - 45 ) 

then in the limit A — > oo the function <r(t) satisfies the a-form of P t , namely, 

(a") 2 = 2a- 2ta' - 4 (a') 3 • (2.46) 
Also, if (2.45) is replaced with 

q m , (t in ,)=-4Af+4A 2 g(t) ; (2.47) 
the limiting equation for q(t) is P l in the conventional form: 

q" = Qq 2 + t. (2.48) 



Analytic properties 

The only branch points of P VI and P , transcendents in the extended complex t-plane 
are 0, 1, oo and 0, oo, respectively. The corresponding tau functions are holomorphic on 
the universal covers of P 1 \{0, 1, oo} and P 1 \{0, oo}. The functions q and o may also have 
movable poles associated to zeros of r. We introduce the branch cuts (— oo,0] U [l,oo) 
(for P ' ) and (— oo, 0] (for P , ) and adopt the principal branch convention for all 
fractional powers of t and 1 — t. 

3. Conformal blocks and instanton partition functions 

3.1. Conformal blocks 



Here we review basic notions about conformal blocks in 2D CFT [5]. For the sake of 
brevity, simplicity and relevance for the rest of the presentation, we will concentrate on 
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conformal blocks for the 4-point correlator on P 1 and Virasoro algebra 



[L 



ai Lm\ 



in 



m) L n+m + — (n 3 



n 



)S. 



n+m,0- 



Only chiral primary fields O will be considered, i.e. their antiholomorphic conformal 
dimensions A^ = 0. We do not require invariance of correlators under the braid group 
action on the positions of fields to avoid constraints on dimensions. 

The three-point correlator of primary fields is fixed by conformal symmetry up to a 
constant factor, 



(0 3 (z 3 )0 2 (z 2 )0 1 (z 1 )) 



rVA A A ^ ~A 3 -Ai-A 2 A1-A2-A3 A2-A1-A3 
^ 1^3, ^2, Ag Z 2 i Z 32 



"31 



where 



Zi — Zj and Aj stand for holomorphic dimensions. Thanks to conformal 
invariance, it suffices to consider more special coordinate dependence, namely, we 
may set Z\ = 0, z 2 = t, Z3 = R with R ^> t. It is also customary to define 
(0(oo)...) = lim i? ^ OOJ R 2A o(C( J R)...), so that, for instance, (O 3 (oo)O 2 (t)Oi(0)> = 
C(A 3 ,A 2 ,A 1 )^- A i-^. 

Besides primary fields, conformal field theory also contains their descendants 
L-\0 = L-\ N . . . L-\ x O, naturally labeled by partitions A = {Ai > A2 > . . . > Aat > 0}. 
Partitions can be identified in the obvious way with Young diagrams. As they play an 
important role in the rest of the paper, we take the opportunity to fix some notation for 
later purposes. The set of all Young diagrams will be denoted by Y. For A G Y, A' denotes 
the transposed diagram, Aj and A' the number of boxes in zth row and jth column of A, 
and |A| the total number of boxes. Given a box G A, its hook length is defined as 
h x (i,j) = \ i + \' j -i-j + l (see Fig. 2). 



/v, 3 J 


|A,|=15 

















































X / ,—X 5 —2 



Fig. 2: Young diagram associated to the partition A = {6, 5, 3, 1} 



Conformal symmetry also allows to compute in explicit form the three-point functions 
involving one descendant: 

(0 3 (00) 2 (t)£-A0i (0)) = C (A 3 , A 2 , Ax) 7A (A 1; A 2 , A 3 ) r™, (3.1) 
(L_ A 3 (00) 2 (t)0i (0)) = C (A 3 , A 2 , Ax) 7A (A 3 , A 2 , Ax) t™, (3.2) 
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where [38] 



7A (A 1; A 2 ,A, 



N 
3=1 



A ? 



A,A 2 



fc=l 



(3.3) 



The action of L„a on the field at infinity in (3.2) should be understood as a result of 
successive contour integration with the energy-momentum tensor. This action can be 
transferred to the fields at and t by deformation of the contour. 

The relation (3.2) is extremely important as it allows to determine the coefficients of 
the operator product expansion (OPE) of primary fields 

2 (t)0i(O) = E C ( A - A2 ' Al ) ^ ( A °> A2 ' Al ) t A ^- A ^L^O a (0) . (3.4) 

Indeed, assuming orthonormality of the basis of primaries, (O a (oo) O/j (0)) = 5 a/ g, 
and considering the correlator of both sides of the last relation with the descendant 
L_\O a (00), one finds that 

(3 X (A a , A 2 , Ai) = J2lQ (Aa)]^ ( A -> A 2, Ax) , (3.5) 

where Qa^ ( A a) = (L-\O a (00) L_ /1 O a (0)) is the Kac-Shapovalov matrix. It can be 
computed algebraically as the matrix element of descendant states 

Q\[i (A) = (a\l Xi ...l Xn l |A), (3.6) 

where |A) and (A| denote the highest weight vectors annihilated by all L n> o and, 
respectively, all L n<0 , satisfying L |A) = A|A), (A|L = (A|A and normalized as 
(A|A) = 1. It is easy to understand that Q(A) has a block-diagonal structure: 
Q\fM (A) ~ 5\\\m- 

We can now finally calculate the four-point correlator (Oi(oo)O 3 (l)O2(t)Oi(0)) . 
Replace therein the product of fields O2(i)Oi(0) by the OPE (3.4) and then use (3.1) 
and (3.5). The result is 

(©4(oo)0 3 (l)02(t)0i(O)) = (3.7) 
= c ( A 4> A 3, A Q ) C (A a , A 2 , Ai) t A *- A ^ A > T c (A 1; A 2 , A 3 , A 4 , A a ; t) , 



where we have introduced the notation 

J- c (A 1 ,A 2 ,A 3 ,A 4 ,A;t)= ^ 7a (A, A 3 , A 4 ) [Q (A)]^7 m (A, A 2 , Ai) t |A| . (3.8) 

A,/xeY 



The representation (3.7) separates model-dependent information (three-point functions 
C (Aj, Aj, Afe)) from the universal pieces fixed solely by Virasoro symmetry. 

The function (3.8) is called four-point conformal block. It is a power series in t with 
coefficients depending on four external dimensions Ai j2j 3 5 4, one intermediate dimension A, 
and the central charge c which enters via the Kac-Shapovalov matrix. These coefficients 
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can in principle be calculated using (3.3) and (3.6). For the reader's convenience, we 



reproduce below several first terms of the series: 



T C (A 1 , A 2 , A 3 ,A 4 ,A;i) = 1 + 



(A - Ai + A 2 ) (A - A 4 + A 3 
2A 



t + 



+ 



+ 



(A - Ai + A 2 ) (A - Ai + A 2 + 1) (A - A 4 + A 3 ) (A - A 4 + A 3 + 1) 



2A (1 + 2A) 

1 + 2A) (A 1 + A 2 + A ( A - 1 \- + 3 2 ^- A2 ) 2 ) (A 4 + A 3 + 



+ 



A(A-l)-3(A 4 -A 3 )^ 
1+2A 



;i-4A) 2 + (c-l)(l + 2A) 



+ 



Direct (i.e. based on (3.8)) computation of conformal block coefficients becomes 
rather complicated at higher levels. An explicit representation for arbitrary level was 
found only recently in a surprisingly different framework. 



3.2. Af = 2 SUSY theories 

The AGT correspondence [2] relates conformal blocks of 2D CFT to Nekrasov functions 
[43j |H]. These functions represent the instanton parts of ei, e 2 -regularized partition 
functions in 4D Af = 2 SUSY quiver gauge theories. 

The simplest case of AGT correspondence deals with SU(2) gauge theory with extra 
Nf = 2N C = 4 fundamental (i.e. transforming in the spin-| representation of the gauge 
group) matter hypermultiplets with masses /ii,...,/i 4 . Parameters of this theory are 
related to those of 4-point conformal block on the sphere by 



C = 



a 3 — a 4 + 



M2 

/i 4 



ol\ — a 2 + 



a 3 + a A — 



2 

e 

2' 



1 + 



6e 2 



ei£2 
a u (e - a v ) 

ei£2 



v 



£\ + £2, 



1,2,3,4. 



The intermediate dimension A of conformal block is expressed via the eigenvalues ±a of 
the vacuum expectation value of scalar field in the gauge multiplet: 
a (e — a) e 



A 



a 



+ a. 



eie 2 2 

The parameter t (anharmonic ratio of four points on the sphere) in the conformal block 
expansion is related to the bare complex coupling constant t uv on the gauge side by 



t = exp 2mruvi 



Ami 6jjv 
+ 



9 2 uv 



2tt 



Partition function in the regularized theory is an integral over a compactified 
moduli space 9Jt of instantons. The integral is given by a sum of explicitly computable 
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contributions coming from fixed points of a torus action on which are labeled by pairs 
of partitions. On the CFT side, this is interpreted as an existence of a geometrically 
distinguished basis of states in the highest weight representations of the Virasoro algebra. 
For more details, generalizations and further references, the reader is referred to 

When all four masses /ii, 2,3,4 — > 00, the fundamental hypermultiplets decouple and we 
get pure gauge theory. Decoupling only some of them yields asymptotically free theories 
with Nf < 2N C . From the gauge theory point of view, the parameter £jy. can be considered 
as a dynamically generated scale. Taking into account the RG dependence of the coupling 
constant, one finds that in the decoupling process this scale should transform in the 
appropriate way: fi^ f — > 00, t Nf — > 0, £jv/-i — l^Nt^Nr fixed [THl E2] • The corresponding 
Nekrasov functions are related to irregular conformal blocks [25| EHl [39] . 

Conformal blocks relevant to Painleve VI equation (27] are characterized by the 
central charge c = 1 and external dimensions 2 [y = 0, t, l,oo) so that we can set 
ei = -e 2 = 1, cti = 9 , a 2 = 9 t , a 3 = 0\, a 4 = 6^ and 

/J, 1 — 9 1 — Ooo, /i 2 = 8 — 8 t , ^3 = 00 + Qt, ^4 = 01+ 0oo- 

Under such identification of parameters, the scaling limits corresponding to the 1st 
line of the coalescence scheme in Fig. 1 describe successive decoupling of the matter 
hypermultiplets: 

Nf = 4 M -+ 00 Nf = 3 -> 00 Nf = 2 M 2 -> 00 Nf = 1 w -> °= pure gauge theory 

OPvi) > (^v) > ' ( P m 2 ) > ( P m 3 ) 

Fig. 3: Decoupling of matter hypermultiplets 

This observation will be used in the next section for the construction of combinatorial 
series for P.. and P,„ tau functions. 

v m l,2,3 

4. Solutions 

4-1- Painleve VI 

Let us first recall the main result of [27J as well as some motivation and evidence for it. 

Conjecture 1. Generic P Yl tau function can be written in the form of conformal 
expansion around the critical point t = 0: 

V W = J2 C ^ fa, 0u 0i, 0oo, ct + n) si & + ^ 2 - e o- 6 'B vl (0 , t , 0^0^, a + n; t) . (4.1) 

The parameters a and s VI play the role of two integration constants, B vl (0q, 0t, 0i,0oo, o~; t) 
coincides with conformal block function T c =\ [0\, 0\, 0\, 0^, a 2 ; t) and is explicitly given by 
combinatorial series 

B m (0o,0 t ,0i,0oo,o-,t) = (l-tf e ^ B^(0 o ,0t,0x,0o O ,cT)t^, (4.2) 
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B 



= n 

>< n 



(7 



#0 2 ) ((*! 



(T 



^(z,j) (A; + ^ - i - j + 1 + 2<r) 
((<9 t - o- + » - j) 2 - gg) ((^i-^ + ^-jf-C) 



(4.3) 



JTie structure constants in (4-1) are given by 

n g \ i + e t + ee - 

e e'=± 

Cyi (^0, #1, #00, 0" ) = ~ 



e'er, 1 



eV 



where G 



(4.4) 



and denotes the Barnes function (see 



ttx, . . . , ct n 
Appendix A ). 

The above claim was obtained in [27J by identifying r VI (i) with a chiral correlator 
(C , £ oo (oo)(9£ 1 (l)(9£ t (t)(9£ (0)) of Virasoro primary fields indexed by matrices L v G 5t 2 (C) 
which are related to monodromy matrices of the auxiliary linear problem for P VI by 
M. v = exp2niC u G SL (2, C). The dimensions of Oc v are equal to A u = |Tr C? v = 9%. 

Our main assumption is that this set of primaries closes under OPE algebra. The 
conservation of monodromy then implies that the dimension spectrum of primary fields 
appearing in the OPE of Oc t (t)Oc (0) is discrete and has the form (a ot + Z,) 2 , where 
2cos27rcrot = Tr M. a M. t . This fixes the structure of the s-channel expansion (4.1) upon 
identification a = a Qt . The constants C VI are obtained from Jimbo's asymptotic formula 



[3T] interpreted as a recursion relation, whereas (4.2)-(4.3) is nothing but the AGT 
representation for c = 1 conformal block rewritten in terms of P VI parameters. 

In fact, Jimbo's formula also expresses the second integration constant s VI in terms 
of monodromy. To give an explicit relation, we need to introduce monodromy invariants 

p v =2 cos 2ix9 v = TiM„, v = 0, l,t, oo, (4.5) 

p fU , = 2cas2n(T lu , = TrMt l M v , fi,v = 0,t,l. (4.6) 
Similarly to the above, the quantities o\t an d cr 01 determine the spectrum of intermediate 
states in the t- and li-channel. The triple a = (aot, o~it, coi) provides the most symmetric 
way to label P VI transcendents with the same 9 = (9 , 9 t , 9i, 6^). The elements of this 
triple are not independent: they satisfy a constraint 

4, (4.7) 



PotPltPO! + Pit + Pu + Pol - ^OtPot - U lt p\ t - ^oiPoi + ^4 



where 



wot = PoPt +PlPoo, 
Ult = PtPl+PoPoo, 
WOl = PoPl+PtPoc, 

uj A = p 2 + p 2 t + p\ + /& + PoPtPiPo 
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Hence, for fixed o"ot, Pit there are at most two possible values for p i- 
Now s VI can be written as 

= (P'u ~ Pit) ~ (gm ~ goi) e2 " CT0t 

VI (2 cos 27r (9 t - o-Qt) - Po) (2 cos 2ti (9 x - a ot ) - Poo) 

where we have introduced the notation 

p'ot = u ot -Pot-PitPoi, 
p'u = -Pu-PotPou 
Poi = Uoi -Poi -PwPit- 



(4- 



Remark 2. Combinatorial expansions of type (4.1) can also be found around the two 



remaining critical points t = 1, oo, as their role is completely analogous to that of t = 0. 
For instance, the series around t = 1 is obtained by the exchange 



This gives 



Xai{0,a)r m {t) 



cot -H- au, 



Poi <->Poi- 



(4.9) 



J2 C Yl (6 h Q u Oo, ^oo, o u + n) S - (1 - t)^ +n ^' 9 ' B, 



(4.10) 

hi Ot, Oo, Ooo, 0\ t + n; 1 — t) , 



with B VI and C VI defined in (4.2)-(4.4) and 



(Pot -Pot) - (p'oi-Poi)e" 



-27riait 



(2 cos 27r (6>i - axt) - pi) (2 cos 2ir (9 - a u ) 



Pc 



(4.11) 



Since the normalization of r vl (t) is already implicitly fixed by (4.1), the expansion (4.10) 

-l 



contains an additional overall constant factor xoi{0,&) . Finding explicit form of this 
connection coefficient is an important open problem which will be treated in a separate 
paper. Note, however, that Xoi(0,cr) disappears from P VI functions o~ Yl {t) and q vl (t). 



Painleve VI equation (2.17) allows to compute the tau function expansions near the 



critical points recursively, order by order, starting from the leading asymptotic terms 
determined by Jimbo's formula. Comparing the result with Conjecture [T] provides the 
most straightforward and convincing test of the latter. Keeping all 9, o arbitrary, we 
have checked in this way (see Section 3 of [27] for the details) about 30 first terms of 
the asymptotic expansion. Also, in a few special cases where P vl solutions are known 
explicitly, the check can be carried out to arbitrary order. This includes Picard elliptic 
solutions and the simplest solutions of Riccati/Chazy type, which correspond to Ashkin- 
Teller conformal blocks [51] and correlators involving low-level degenerate fields. More 



complicated Riccati solutions are discussed in Subsection 4.3 of the present paper 



Numerical efficiency of the expansions (4.1 ), (4.10 ) is illustrated in Fig. 4. For random 



complex 9, a we plot on the same graph the series for o~ vl (t) around t = (blue line) and 
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Re(ovi) 




Re(o-vi) 



L= 15 




Im(o-vi) 




Im(o-vi) 



L= 15 




Fig. 4: Truncated P VI series at t = and t — 1 with 
(0.1902 + 0.3106i, 0.4182 - 0.2109i, 0.3429 + 0.3314i, 0.0163 + 0.1805i), 
a = (-0.3272 - 0.481H, 0.0958 + 0.3168i, 0.4762 + 0.1028z) 



t = I (red line) keeping the terms up to O (t L ) and O y(l — t)j with L = 5, 15. Zooming 
near the endpoints t — 0, 1 would display oscillations of rapidly increasing frequency and 
decreasing amplitude due to non-zero imaginary parts of c%, a u . 



4-2. Painleve V and Ill's 



Next we consider the scaling limit P VI — » P v given by ( 2.30 )-( 2.33). Conformal 

cr; I") has a well-defined limit as A — > oo, which 



block function B VI (6 Q ,6 t , ^ 

can be calculated termwise in (4.2)-(|4.3|). The asymptotics of the structure constants 
A+e* a- 



— , cr) ensures consistency of the expansion (4.1) with the limit ( |2.33 ) 

;h 





vi V u ' £ ' 2 ' 2 

More precisely, using the estimate (A.l ) from the Appendix A, it is easy to check that 
lim ^ _<T ^ VI ^' ^2^' ~ 

A— >oo 



G 2 (1 + A) 



C v (^O) #t) °") 



(4.12) 



where 



C v (^0) #t ; ^*) o") 



IM 1 



ecr, 1 



e=± 



?o + ecr, 1 
-2ea 



0t - 0n + ecr 



(4.13) 
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One could even completely get rid of the denominator in the l.h.s. of (4.12) by modifying 
the normalization of P VI tau function in (4.1) (e.g. by dividing all structure constants in 
(4.4) by a cr-independent factor G 2 (1 + 9\ + 6^)). 
Altogether, this leads to 

Conjecture 3. P v is solved by the following tau function expansion at t = 0: 

r v (t) = J2°v (0o, e t , 6*,o- + n)s n v & +n ? B v (9 , 9 t , 9*,a + n;t). (4.14) 

Here again a and s v are arbitrary parameters, irregular conformal block B v (9q, 9 t , 9*, a; t) 
is a power series defined by 

B v (9 ,9 t ,9*,a;t) = e- e < t ]T (0 O A * |A|+ImI , (4.15) 



B^(9 ,9 t ,9*,a) 



+ a + i-j)((9 t + o- + t-j) 2 -9 2 ) 



n 



x (4.16) 



/* — cr 



X 



j) ((* 



t ~ O 



3 



(J)L hl{i,j)(X i + ti , j -i-j + l-2<ry 



01) 



and the structure constants C v (9o,9t,9*,o~) are given by (4- IS). 

The second P v critical point t = oo corresponds to irregular singularity of the 
associated 2x2 linear system and is obtained by the fusion of two P vl critical points 
1, oo. The expansion around this point cannot be extracted from P VI series and requires 
the knowledge of complete irregular OPEs. For the same reason, we are so far unable 
to treat P lv , P n and P r However, long-distance expansions of this kind are available in 
a few special cases where the solutions of P vul can be expressed in terms of Fredholm 
determinants, see Section [5] 

Because of the presence of irregular singular points, monodromy data for P v involve 
Stokes multipliers. The expression for the integration constants a, s v of Conjecture [3] in 
terms of monodromy can be extracted from Jimbo's paper |31j . 



Repeating the previous arguments almost literally for the scaling limits ( 2.34 )-( 2.37), 



(2.38)-(2.40) and (2.41)-(2.43), one obtains short-distance expansions for tau functions 
of three nontrivial P m equations: 

Conjecture 4. Expansion of r n , (t) at t = can be written as 
where the irregular conformal block B , (9*, 9*, o~;t) is given by 

n _w|A|+H 



(4.17) 



B m , (0*A, a;t) = c ~i ]T B^ 



(4.18) 
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6* + cr + i-j)(6* + o- + i-j) 



x 



n 
n 



(i)ex hl&fi {X' j + f H-i-j + l + 2a)' 
9* — a + i — j) (9* — a + i — j) 



x 



and the structure constants can be written as 

l + 9* + ea,l + 9+ + ea 
l + 2ea 



2 ■ 



c m , (e„e„a) 



n«[ 



Conjecture 5. Expansion of T m , (t) at t = is given by 

^ (*) = E C n:^ * + «) C/^X* * + "5 *) 

arbitrary a, s m , and 

B m ,{9^a-t)= bQ> (0.,*)tM+M 



= n 



o- + z - J 



X 



9* - a + i - j 



n 



^ (^^)=n 



G(l + 9, + ea) 



2 ■ 



G(l + 2ea) 



Conjecture 6. Expansion of P m , tau function at t = is: 

T (t) = \"C ,(a + n)s n & +n)2 B ,(a + n:t) 

in' V / / j in' v / ni; hi' v ' 



where 



B 



(in' 3 ) 
A,/t 



n -2 



16 
(4.19) 



(4.20) 
(4.21) 

(4.22) 
(4.23) 



(4.24) 

(4.25) 

(4.26) 
(4.27) 



» = [ II ^frj) (A; + ^ - i - j + 1 + 2a) 

(m)£A 

x II h fl (i,j)(\ i + f i' J -i-j + l-2o-) 

C m ,(a) = [G(l+2a)G(l-2a)}-\ (4.28) 

As in the P m case, Conjectures [3jj6] can be verified by iterative reconstruction of the tau 
function expansions from the leading asymptotic terms using the equations P , 
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4-3. Classical solutions: AGT vs Gessel's theorem 



In 2002, Forrester and Witte [24] have proved a remarkable determinant representation 
for a family of Riccati solutions of P VI . Their result can be restated as follows. Define a 
five-parameter family oi N x N Toeplitz determinants 

N-l 



D 



N 



(t) = det 



J j,k=0 



(4.29) 



[1-ty 



T (1 + T] — m) r (1 — rj + m + v' 



r (1 — rj + m) r (1 + rj 

Then the function 



m 



v) 



■ 2F1 
2F1 





t 1 


. 1 + 7] — m 


t - lJ 


- 1 + v, -v' 


t - 


- 1 — i] + m 


t - lJ 



+ 



iV(iV+i/+i/) 



is a tau function of P VI with parameters 

1 



D 



N 



(t) 



(4.30) 



(4.31) 



?0, ^t, ^1, 6*c 



(77, AT, —N - v 



v' + rj) . 



Looking at the asymptotic expansions oiD%' u '^\t) atOandl, one can also identify 



the monodromy exponents 
(a t, Cit, CTqi) 



(AT + 7], v + v\ N + v - v + r]) 



Almost all structure constants in (4.1 ) vanish because of the relations 9t — %, Uot — do + 9t 
(recall that Barnes G-function has zeros at negative integer values of the argument). The 
only non-zero constants correspond to n = 0, —1, . . . , — N, so that there remain only A^ + l 
conformal blocks. The parameter s VI in (4.1) is related to £ in (4.30) by 

sin7rz/sin7r(r] — v') 



sin ttv' sin 7r(r] + v) 



(4.32) 



Let us now consider in more detail the case £ — > 0. Then (4.32) implies that s VI — > 00, 
which means that the expansion (4.1) at t = contains only one [n = 0) conformal block 
E m (f ' f > E± T^ L , ^ ±R , ^y 2 ; tj . The product over boxes of fi in the AGT representation 
(4.3) contains a factor i — j due to the relation a ot = 9 + 9 t . Since this expression 
vanishes for the box (1,1), the quantity does so for any non-empty /i. Moreover, 
the factor i — j + N in the product over boxes of A reduces the summation in (4.2) to 
Young diagrams with Ai < A" (i.e. with the length of their first row not exceeding N). 
Therefore, Conjecture [T] for the above parameters is equivalent to the following identity: 

i-j + N (i-j-v)(i-j-i/ + rj) 



D 



N 



c N J2 1 

AeY|Ai<iV 



|A|- 



Nri 
2 



n 



j + N + r] 



h 2 



(4.33) 
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(4.34) 



where the constant prefactor 

r 1 + N, 1 + v' + N, 1 + T], 1 - 7] + u' 
N ~ l 1 + t] + N, 1 - 77 + v' + N, 1 + v' 

can be computed using pure Fisher-Hartwig determinant. 

In the limit 77 — > 0, the left hand side of (4.33) reduces to N x N Toeplitz determinant 
with the symbol 



a(o = (i+vicy (l+vtc'Y . 



(4.35) 



Also, Cjy = 1 and the first factor in the product on the right disappears so that the 
r.h.s. coincides with the length distribution function of the first row of a random Young 
diagram distributed according to the so-called ^-measure [8]. The equality (4.33) can then 
be rigorously demonstrated using a dual version of Gessel's theorem [221 EZ]- 

Remark 7. We draw the reader's attention to the fact that Toeplitz determinant with the 
symbol (4.35) with v = —v' = | coincides with diagonal two-point Ising spin correlation 
function on the infinite square lattice. Its relation to P VI is rather well-known [32]. It is 
intriguing, however, that this lattice correlator is equal to a (particular limit of) conformal 
block in continuous 2D CFT with c = 1. 

Analogous results for P v and P m , can be obtained by successively sending v' and v 
to infinity. For instance, consider instead of A^ u ' v ' ,v '^ (t) and t^' v ' v '^(t) the quantities 



t~ 



Fx (-is, 1 + 77 



T (1 + i] — m 

^(l + ^t^V* 
T (1 — 77 + m) r (1 + 77 — m + u) 



m 



-t) + 
xFi (1 + v, 1 



(4.36) 

77 + m,t) , 

(4.37) 

then t^' (t) is a tau function of P v with (6 , 6 t , #*) v = \ (77, N, N + rj + 1v\ Similarly, 
if we define 



T p v '®(t) =t^^e^ det 



.4 



j-k 



-I N-l 



j,k=0 



L r)—m 



(2Vt) + £I TO _„ (2^) , 



T 



N 



(t) = t~ 



"2 det 



then Ttf'®(t) is a P , tau function with 



N-l 

j,k=0 ' 

N+ri n 
2 ' a * 



(4.38) 
(4.39) 



N-r) 



For £ = and 77 







the symbols of Toeplitz determinants (4.37), (4.39) are smooth and can be written as 
(1 + ^tCYe^- 1 and e^( f+c_1 ). Gessel representations of these determinants coincide 
with the results derived from Conjectures [3] and [4} 

In the general case £ ^ 0, the function T (^',v,0(t) 

is a polynomial of degree N 

in £. The coefficients of iV + 1 different powers of £ are s-channel conformal blocks 
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with internal dimensions (6q + @t — k) 2 , where k — 0, . . . , N. Alternatively, one can first 
transform hypergeometric functions to make them depend on 1 — t and then expand the 
determinant in powers of £, the analog of parameter £. The result has the form (4.10) 
of a sum of t-channel conformal blocks with internal dimensions (6\ + 6 t — k) 2 , again 
with k — 0, . . . , N. The relations between the expansion parameters are given by (4.32), 

iKi = t s vi = £ and ( x - ^viX 1 - Ki) = 1 + £• 

The CFT interpretation of this picture is as follows. The tau function (4.31) is 
a four-point correlator of primaries which involves level N + 1 degenerate field (here 
Oc t (t)). Its expansions at t = and t = 1 incorporate all allowed intermediate dimensions. 
Determinant representation (4.29) can in fact be used to compute the fusion matrix for 
the corresponding two sets of conformal blocks. This task simplifies in the case £ = 0, 
where we are left with one s-channel block tranforming into a linear combination of iV + 1 
t-channel ones. 



5. Examples and applications 

5.1. Integrable kernels 

In many applications of Painleve equations the relevant tau functions can be written as 
Fredholm determinants of scalar integral operators of the form det (1 — K\i), where K\i 
denotes the restriction of the kernel K(x, y) to some interval I C R. These kernels usually 
have integrable form, that is 

K(x,y) = X ^ {y) -^ xMy \ A G C. (5.1) 
x-y 

As is well-known, given / = (Jy=i i a 2j~i, a 2j) and (p, ip verifying the differentiation 
formulas 

( I'd ) = A{x) ( lU ) ' 

with some rational matrix A(x), the corresponding Fredholm determinant satisfies a 
system of PDEs with respect to {aj} [56J. For <p, %p given by classical special functions 
and sufficiently simple /, this system can often be solved in terms of Painleve functions 
[231 ESI EH ESI EHl [59]. 



5.1.1. Hypergeometric kernel. The most general known example corresponds to the 
choice 



<Pg( x ) 



1 + u + rj, 1 + v' + r( 
2 + v + v' + 77 + 77' J (]_ 



X 2 



X) 2 



2^1 



1 + v + 77', 1 + v' + 77' 

2 + Z/ + V' + 7] + 7]' 



X 



X — 1 - 
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1 + v + rf, 1 + v' + 7] 
l + v + v' + rj + r)' 



X 2 

. v+v' + 2ri' 
. — X) 2 



v + 77', i/' + 77' 

+ v' + r/ + 7]' 



X 



x — 1 



with A = 7r 2 sin7rz/sin7rz/. The kernel K Q (x,y) contains four parameters z/,z/, 77, 77' G C 
chosen so that the Fredholm determinant 

D G (t) = det(l-K G \ m ), f e(0,l). (5.2) 

is well-defined. We will not try to determine the set of all possible values of v,v f ,77,7/; 
the interested reader may find examples of admissible domains in |10j . 

The above 2 -Pi kernel first appeared in the harmonic analysis on the infinite- 



dimensional unitary group pm Ej- Later it was shown |35j that the determinant (5.2) 
coincides with a correlator of twist fields in the massive Dirac theory on the hyperbolic 
disk [T71 EU HE] . From the point of view of the present paper, the most interesting feature 
of D G (t) is that it is a Painleve VI tau function, see [10] and also Sec. 5 of [35] for a 
simpler proof. The corresponding P VI parameters are 

(6» , 9 t , 0i, 6>oo) vi = - (u + v' + 7] + r)', 0, v - v\ rj - rf) . 

Monodromy characterizing this particular solution is determined by 
v + v' + 77 + 7/ v + v' 



cot 



fit 



2 ' " 2 ' 

cos27rcroi = 2e _7ri< - r?+r? ) sin7rz/sin7rz/' + cos7r(r7 — 7]\ 



Substituting these parameters into (4.11), it can be easily checked that s 



= 1. 

Remark [2] then implies that the large gap (t — > 1) expansion of the 2 -Pi kernel determinant 
is given by 

(5.3) 



D g(*) = ^ C G (z/ + n, z/' + 72, 77 - n, 77' - n) (1 - t) 
x B G (v + n, v' + n, 7] — n, 7]' — n; 1 — t) , 



where 



C G (z/, z/, 77, 77') = G [1 + 77, 1 + 77', 1 + 77 + v + z/, 1 + 77' + v + z/] Yl G [ 



e=± 



1 + ez/, 1 + ez/' 
+ e (z/ + z/) 



r , / / , ,, r ,'v-v' + + 77 + V 77-77' Z/ + Z/' 
B G (z/, z/ , 77, 77 ; 1 - t) = B, r , ( — - — , 0, , — - — , — - — ; 1-t 



and £> VI is given by (4.2)-(4.3). Also, [351 Conjecture 8] suggests that the constant x G is 
equal to 



X G = G 1 + 77 + z/, 1 + 77 + z/, 1 + 77' + z/, 1 + 77' + z/ 



(5.4) 



Constructing the expansion at t — is less straightforward. It is of course possible 
to compute a few first terms in the small gap asymptotics directly by expanding D G (t) 
into Fredholm series. This yields, for instance, 

D G {t) = l-K G t l+ ^' + ^' [l + o(l)], 
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a r 



with 

1 + 77 + V, 1 + T]' + V, 1 + T] + is', 1 + 77' + V 1 

2 + ri + ri / + v + v',2 + ri + ri' + v + v' 

On the other hand, direct application of Conjecture [T] is ambiguous because of special 
parameter values. First, Barnes functions G (1 + 6 t ± (0 O — Pot — n )) in the structure 
constants vanish for n ^ 0. At the same time the quantity s VI diverges due to zero 
denominator. The right way to handle this is to fix the values of 0's and a±t, and then 



consider the limit o~ot — > #0 with the help of the formulas (A. 2 )-( A. 3 ) from the Appendix A. 



C G (v, u',r],r}',n) = (-A) n G 



x G 



(5.5) 



The result is that only the terms with n > survive in the sum over n and the structure 
constants reduce to 

l + nA+n + n' + v + v' + n ~| 2 

x 

1 + r] + T}' + v + u' + 2n 

I + 77 + Z/ + n, 1 + 7]' + v + n, 1 + 77 + z/ + n, I + 7/ + v' + n 
1 + 77 + v, 1 + T]' + v, 1 + 77 + is', 1 + 77' + v' 

In addition, because of the factors 9 t + i — j ± (er t + n — 9 ) in the products over boxes 
of A, n G Y combinatorial summation in conformal blocks can be restricted to Young 
diagrams with Ai < n, < n. 

This leads to the following expansion of D G (t) near t = 0: 



(5.6) 



n=0 



where C G (u, u', 77, 77', n) is given by (5.5) and 



=n 



(i — j + n) (i — j + n + 77 + 77 + v 



j + n + 77 + v) {% — j + n + 77' + v) 



h 2 x {i,j) (A£ + ^ - i - j + 1 + 2n + 77 + 77' + v + v')' 



x 



j — n) {1 — j — n — 77 — 77 



v — v 



n 



v 



n 



T]' - V' 



h l(hj) (Xi + fJ.j-i-j + l-2n-r]-rf 



Note that individual conformal blocks in the sum over n in (5.6) give the corresponding 



terms in the Fredholm series of the 2 -Pi kernel determinant. Numerical checks for randomly 



chosen 77, 77', u, u' show that the expansions (5.3 ) and (5.6 ) perfectly match for intermediate 



values of t. In particular, this confirms the conjectural expression (5.4) 



5.1.2. Whittaker kernel. The Whittaker kernel P, [TO] emerges in the limit 



K w {x,y) 



lim -K c (l--,l-^- 
77' V 77' 77 
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It contains three parameters u,u',r} and has integrable form (5.1 ), A is the same as above 
and 

!f w (x)= r (1 + 77 + V) X~^W v+ u'+2r, . 1 v -u' (X) , 

2 ' 2 ' 2 

tfj w (x) = I (1 + T] + U') X~*W v+v > + 2r, 1 (x) , 

2 2 ' 2 

where W4,m (x) denote the Whittaker functions. 
Fredholm determinant 

D w (t) = det (1 -K w \ {t)0o) ) , te{0,oo), 

is related to a particular Painleve V tau function by 



D w (t)=t —T v (t), 

Its expansion around t = may be found from 

t 



(5.7) 
(5.8) 



D w (t) 



lim L> 1 

r/'^oo e ' V 77' 



Namely, the appropriate termwise limit of (5.3) gives 



(5.9) 



where the limits of structure constants and conformal blocks are 

1 + i/,1-i/,1 + z/,1-z/,1 + J7,1 + »7 + i/ + i/ 
1 + v + u', 1 - v - z/, 1 + 7] + v, 1 + rj + v' 



v — v 1 v + v' v + v' 



and fi v was defined in (4.15)-(4.16). Although we are not able to write similar 



combinatorial expansion at t = oo, in the latter case D w (t) can still be expanded into 
Fredholm series. Hence, for example, 

D w (t -> oo) = 1 - A T (1 + 77 + v) V (1 + rj + v') e -t r (2+2 V+ ,+,') ^ + Q ^-i^ 

5.1.3. Confluent hypergeometric kernel. Another interesting scaling limit of the 2F1 
kernel corresponds to setting 

v' — i/ — iA, V = Vo + 
and then considering 



I / x 

KJx,y) = lim — K c I — , 
FV ' a^oo A G VA 



x y_ 
A 
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The result is the so-called confluent hypergeometric kernel JTOJ, [15] . It depends on three 
parameters 



r+ = u + rj', r_ = u' + rjo, £ 



1-e 



2niu 



T ( r __ r+ ) 



27T 



and is given by (5.1) with 

1 + r+, 1 + r_ 
1 + r + + r_ , 2 + r + + r 



^ p (x) = X' 



e 2 i-Fi (r + + l,r+ + r_ + 2, ix) 



ip v (x) = x 2 e 2 i-Fi (r + , r + + r_, ix) . 



Similarly to (5.7)-(5.8), the i_Fi kernel determinant 



D F (t) = det(l-ir p | (Q!t) ), te(0,oo), 
can be expressed [10] in terms of a Painleve V tau function: 

(r + +r_) 2 

D F (t)=t *— r v (it), 



1 



(r+ + r_, 0, r + — r_" 



(5.10) 
(5.11) 



Note that D F (t) = limA-s.oo D G [J^ J ■ Applying this termwise to (5.6) and using the 
properties ( |A.1[ )-(A.2) of the Barnes function, we derive the expansion of D F (t) at t = 0: 



n=0 A,/ieY|Ai,Ati<n 

where 



C F (r+,r_,n) = G 



1 + n, 1 + r+ + r_ + n 
1 + r + + r_ + 2n 



n 2 



G 



1 + r + + n, 1 + r_ + n 
1 + r + , 1 + r_ 



t-t (i — j + n) (i — j + n + r+) (i — j + n + r + + r_) 



n 



(i i)6 A ^a(m) W + A** - * - J + 1 + 2n + r + + r_)' 
(z — j — n) (i — j — n — r_) (i — j — n — r+ — r_) 
^(m) + ^ - « - J + 1 - 2n - r+ - r_) 2 



x 



5.1.4- Sine kernel. Certain specializations of the i_Fi kernel play an important role in 
random matrix theory. In particular, for r + = r_ = r it coincides with the Bessel kernel 

[SB E2| 

r£v^+I (|) Jr-\ (I) - Jr-l (|) Jr+i (l) 



7T4 



2 x — y 

which in the case r = reduces to the celebrated sine kernel 



K sine (x,y) 



2£sin^ 
x-y 
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It is well-known that the determinant 

-Dsino it) = det (l - # S me | (0,t)) 

for £ = 5- coincides with the scaled gap probability in the bulk of the Gaussian Unitary 
Ensemble [23]. The expansion (5.12) thus gives a complete series for this quantity: 



A,/ieY|Ai,jt4<n 



n=0 



(5.13) 



where 



+ n) 



JJa (A;. + /i,-^-j + l + 2n) : 



x 



n 



[i-j-n) 



First terms of the series (5.13) are recorded in the Appendix B. In particular, they 
reproduce the results obtained by an iterative expansion of the corresponding Painleve V 
solution, cf Eq. (8.114) in [23]. Note that our t = 2iit- 23 -, £ = We have also checked 



the agreement of (5.13) with the known large gap (t — > oo) asymptotics [18 

1 



lysine 



£=JL 



71 G l 



t 4 e 2 



i + -r 2 + — r 4 + o(r 6 ) 

32 2048 V ; 



5. i.5. Modified Bessel kernel. One may also study a further scaling limit of the i_Fi 
kernel by setting 

,7rA 



r re 
r± = - =F «A, £ = £ 



and defining 



K B2 (x,y) = lim -i^ ( 

A— >oo A V 



;2 2vr 

x y_ 
A' A 



Asymptotic properties of the confluent hypergeometric function imply that 

I r+1 (2v^) J r -i (2Vy) " ir-x (2v^) h+i (2y/y) 



x-y 



Fredholm determinant D B2 (t) = det (l — ^ B2 |(o,t)) is related to a tau function of 
Painleve III^ with 9* = —9+ — \ by 

D B2 (t)=t-^eh m ,(t). (5.14) 
Its small gap expansion can be calculated using that D B2 (t) = limA^oc D F m . We find 



oo 



n=0 



1 + n, 1 + r + ra 
r + 2n 



A j( ueY|Ai,/4<n 



(5.15) 
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? B2 i \ TT (i — j + n) (i — j + n + r) 



®zm = n 
n 



d fia h l('-3) ( X 'j + W - i ~ 3 + 1 + 2 " + r ) 
(z — j — n) (i — j — n — r) 



x 



5.1?. Sine- Gordon exponential fields 

A well-known example of appearance of Painleve transcendents in integrable QFT 
is provided by the two-point correlation function of exponential fields Q (mr) = 
(O u (0) O v i{r)) in the sine-Gordon model at the free-fermion point [3151]. The spectrum 
of this model consists of fermionic excitations of mass m, parameterized by the topological 
charge e = ±1 and rapidity 9 6 R. Lattice counterparts of the exponential fields have 
been introduced and studied in [2SJ H7] . 

5.2.1. From form factors to Macdonald kernel. Under normalization (O v ) = 1, the 
exponential fields are completely determined by their two-particle form factors [50] 

F„(9,9 , ) = + -(9;9'\O l/ (0)\vac) 



2tt cosh £=i 



Multiparticle form factors can be written as determinants of two-particle ones. This allows 
to sum up the form factor expansion 

°° 1 POO POO 

Q ( mr ) = — ■■ d9 x ...d9 n e~ mrY ^^ coshefc x 

n=0 ei,...,e n =± ^ ~°° 

x (vac|a(O)|0i, • • -,9 n ) ei j : ^ (6 U . . . ,9 n pA$)\vac) 

to Fredholm determinant Q (mr) = det (1 — K SG ). The corresponding kernel acts on 
L 2 (IR) and is expressed in terms of dressed two-particle form factors: 

/oo 
T- v (9", 9) T v , (9", 9') e - I f : (coshe+2coshe''+coshe')^. 
-oo 

Let us show that K SG (9, 9') is equivalent to a more familiar classical integrable kernel 



K M (x,y) on L 2 ( m ^—,oo). The latter is defined by (5.1) with A = 7r 2 sin7ri>sin7rz/ and 
(p, ip given by Macdonald functions 

<p M (x)=2y/iK v ,-. v+1 (2y/Z), lj> M {x)=2K v ,- v (2yft). (5.16) 

This kernel can be seen as a further scaling limit of the Whittaker kernel from the previous 
subsection. Indeed, one may check that 

K m( x ,V) = lim ~ K w ( ~> - 

By equivalence of K SG and K M we mean that Ty K™ = TrK™ for any n G Z> . 
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Here is a proof. First note that the Macdonald kernel admits an alternative simple 



K M ( x ^y) = x / ^ M ( xt )^M (y f ) dt - 



form 



i 

This representation results from the identity 
d 



(5.17) 



which is itself an easy consequence of the differentiation formulas 



x 



dx I ip M (x) 



x \ ( Vmi x ) 



On the other hand, parameterizing the rapidities as« = e e , one can write K n = Tr K™ as 



** . . . || || (»■>+»;')} , (5 . 18) 



j=l i= l 



Uj + Uj+i 



with «2n+i = U\. Now make in (5.18) the following replacements: 



e 2 



'(U2j_l+«2j) 



- t v- l( - u v- 1+U2 ^dt 2j -i, 



i n r 

2 



mr 

e 2 



(m 2 /+m 2 / +1 ) 



U 2 j-l + U 2 j 



U 2j U 2j+l 



-t 2j {u 2j +u 2j+1 ) dt 



2 



where j = 1, . . . , n. This yields a 4n-fold integral 



dui 
o u i 



du 



2n 



U2r, 



dti 



m r 
2 



m r 
2 



2j-l, 



- 1 2 j i- 2 U 2 j_ l - 1 2 j i_ i U 2] - 1 2 j U 2 1 



'2j 



i=l U 2j 



with t = ^2n- The variables ui, . . . ,u 2n are now decoupled. Integrating them out with 
the help of the standard integral representation of the Macdonald function 



we finally obtain 



poo roo 2n 

/ dh... dt^TJ^^tj-xtj). 

J mr mr 



«n = A* 

/ mr J mr 

2 2 j=l 

After the change of variables £2,3-1 ^ ^hj-i, ^2j | — ► ~^;*2j the last expression can 



obviously be written as TrK™ with K M given by (5.17). 
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5.2.2. Painleve III and asymptotics. Painleve representations of the two-point function 
of exponential fields [TJ [51] can now be rederived by applying the standard random 
matrix theory techniques [56J to the Macdonald kernel. The final result is that 

coincides, up to a simple prefactor, with a tau function 



det 1 - K 



Q(mr) 

of Painleve 1 11^ equation with parameters 6* = —6+ = ^f-: 

Q(2Vt)=t- (J ^elr n ,(t). 

The integration constants specifying this tau function are [31] 
v + v' 



a 



1. 



(5.19) 



In general, the tau function is defined up to multiplication by a constant, which 
in the case at hand is fixed by normalization of the VEVs: Q(mr) ~ 1 as r — > oo. 
Subleading corrections to this long-distance behaviour can be obtained from the form 
factor expansion. For instance, taking into account the contribution of two-particle states, 
we find 



1 _ Q(2Vt) = A r(v' M (xW M (x) - <p u {xW u {xj) dx+O (t 



-1 _-8%/t 



Short-distance asymptotics of Q(mr) is also known. Assume that |Re {y + u') | < 1, then, 



as t — > 0, 



Q(2v / t) 



The value of a in (5.19) is determined by the exponent vv\ found in [31J. The coefficient 
C SG (u, v') was calculated by Basor and Tracy in [3]: 

C 5G {vy) = G\ l + » A ~ U \\ +U ' A ~, U '\ (5-20) 
SG v ' L 1 + v + u', 1 - v - v' J v ' 

Note that the last expression coincides with C m , (^f-, ^y^, ^2^) defined by (4.20). This 
simply means that the normalization of r , (t) in Conjecture 4 corresponds to setting 
(O v ) = 1 in the sine-Gordon case. 



5.2.3. Short- distance expansion of {O u (0) O v i{r)). We are now ready to write complete 
short-distance expansion of the two-point correlator Q(mr). Combining the above with 
Conjecture [4] gives the following series: 

/ 2 2\ {f+n)(u'+n)+\X\+\n\ 

Q{mr ) = Y^C SG (u + n,u' + n) ^ JB^, {v + n,u' + n) ( ) , (5.21) 
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where C sg (za // 



') is defined by ( |5.20 ) and 



X 



n 

(i,j)ex "a 

n 



(z - j + v) 



j + j/ 



(y + iH-i-j + l + v + iSy 
(i-j - u) (i-j- v') 



h?Ji,j) (Xi + fJj-i-j + l 



>')' 



The series (5.21 ) has a familiar structure of conformal perturbation expansion [191 [62]. The 
non-analytic factors m 2 ( v + n )( v '+ n ) correspond to non-perturbative VEVs of the primary 
fields which appear in the operator product expansion O v (0) G v i{r). All other corrections, 
including the VEVs of descendant fields and CPT, come in integer powers of the 
coupling m 2 . 



Fig. 5 illustrates how well the series (5.21) fits form factor expansion to give all- 
distance behaviour of the correlator. We fix v = 0.3, v' = 0.45 and compute the 
expansion Q 15 (mr) taking into account the terms with n = —4, ... ,4 up to descendant 
level 15, as we did before for P VI . Plots (A), (B), (C), (D) correspond to the logarithms 
of — 2~^=i "Tr-K^Q — hi Q 15 (mr) (solid lines) and jTrK* Q (dotted lines) for £ = 1, 2, 3, 4. 
Hence Q 15 (mr) correctly accounts for the 2-particle form factor contribution to long- 
distance asymptotics up to mr 3.2, 4-particle contribution up to mr w 2.6, 6-particle 
and 8-particle ones up to mr w 2.1 and mr w 1.7. 



-10 



-20 




1.0 



2.0 



2.5 



3.5 



Fig. 5: Short-distance vs form factor expansion 



Let us stress that we are dealing here with a correlation function of massive theory. 
Its description by holomorphic conformal blocks therefore looks rather surprising and is 
presumably related to the affine sl(2) symmetry of the free-fermion sine-Gordon theory 
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5.2.4- Painleve III and 2D polymers. The change of variables q ui , (t) = ^expip(r) 
r = A\fi maps P ' with 9+ = 0, 6* = | to radial sinh-Gordon equation 



r 



sinh2-^. 



(5.22) 



A particular solution of this equation describes universal scaling functions of 2D polymers 
[20| [63] . It is characterized by the boundary conditions 



ib(r -+ 0) ~ -- lnr - - In ^ + O (r 4/3 ) 
^ V ' 3 2 4 V 7 



r 2 (i/3) 
r 2 (2/3) 



6' S m' 



(5.23) 



1. In fact the 



which in our notation correspond to integration constants o~ m/ 

relevant P , solution is a Backlund transform of a solution associated to the tau function 

m i 

considered in the previous subsection. The precise relation between the two quantities is 



sinh 2 ip(r) 



(lnQ(r))" + T- 1 (In Q(r))' 



On the other hand, it is known [16] that P , with 9* = 0, 9* = \ (and hence the 
radial sinh-Gordon equation!) is equivalent to P ' . Namely, if we set 



t 



Qm> 3 [tin's 



then q(t) satisfies appropriate P , . This allows to give an alternative characterization of 

. . m i . . . . 

the solution (5.23) via the expansion (4.25 )— (4.28|): 



-2^(0 



-4r" 



i d d 
— r — 
dr dr 



lnr , (2" 12 r 4 ) 
in' v / 



=1,(7= 



or, in yet another form, 

t- , (2 

III', v 



-12 r 4 l 



=1,(7= 



=1,(7= 



where the normalization of both tau functions in the last formula is precisely the same as 
in Conjecture |6j 

6. Discussion 



We believe that by explaining the title of this paper we have partially answered P. Deift's 
question from the Introduction. Besides the obvious need for rigorous proofs of our claims 
in Section [4], many other questions beg to be addressed. Why instantons? Is there a way 
to obtain irregular "form factor" expansions at oo for general solutions of P v and P in ? 
What about P IV , P u and PJ 
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A particularly interesting problem, already mentioned above, concerns the 
computation of connection coefficients of Painleve tau functions (akin to Dyson-Widom 
constants in random matrix theory). In the P vl case this is very much related to 
determining the fusion matrix for c = 1 generic conformal blocks. 

Another intriguing issue is the quantization of Painleve equations [HJ EOj- The 
existing paradigm usually associates isomonodromic deformations to semiclassical limit 
of CFT [53J. For instance, the scalar Lax pairs for P VI _ : emerge in the c — > oo limit of 
two BPZ-type differential operators [32] ■ The results presented here and in [27] suggest 
a completely different, c = 1 point of view. In this picture, classical Riccati solutions of 
Painleve equations may be naturally deformed to Coulomb /3-integrals. It would be nice 
to understand whether the general case allows for a similar /3-deformation. 
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Appendix A. Barnes function 



Barnes G-function satisfies the functional equation G (1 + z) = T (z) G (z) and is defined 
as the infinite product 



G(l + z) = (2vr)5 e xp 



Z + Z 2 (l+-f) 



n(i+ 



fc=i 



k 



exp 



2k 



where 7 is the Euler's constant, or via the integral representation 



G(l + z) = (2vr)5exp 



dt 

7 



-zt 



\ — e : : " 

4sinh 2 | ~ t + ~2 6 



Rez > — 1. 



It is analytic in the whole complex plane and has the following asymptotic expansion as 
\z\ — > 00, arg z ^ 7r: 

/r 2 1 \ 3r 2 

lnG(l + z) : 



lnz - ^ + Z - hi27r + ('(-1) + O (z~ 2 ) 



~2 ~ 12 

One of the consequences of this asymptotic behaviour is the formula 

1 + z + a, 1 + z — a 
L l + z + P,l + z-P 



G 



z 



a 2 -0 2 



[l + 0(z- 2 )] 



(A.l) 
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Another useful relation is 

1 + z + n, 1 — z 



G 



1 — z — n, 1 + z 



(^— ) , n E Z. (A.2) 



sin TTZ 



It is easy to deduce from it that, as s — > 



n(n — 1) 

G(l+e-ra) ~£ n (-l)^^G(l+n), nGZ> . (A.3) 



Appendix B. Sine kernel conformal blocks 

Consider the functions 



A,/ieY|Ai,/4<n 



|A| + H 



which appear in the expansion (5.13) of the GUE gap probability. Below we record the 



terms contributing to D sinc (£) as at least t 30 : 
£ ino (0;t) = £ inc (l;t) = l, 

» ,„ x t 2 t 4 t 6 t 8 t 10 t 12 

B Binc (2;t) = 1 - — + — - — — + — — — - ^nmooronn + 



75 7840 1134000 219542400 55091836800 17435658240000 

t 14 t 16 t 18 



6802522062336000 3210079038566400000 1803084500809912320000 

^22 

+ 1189192769988708925440000 _ 910206422681575219200000000 + 

j.24 + 26 

+ - - + o (t 28 ) 

800331904605748883816448000000 801284680682660489630515200000000 v ; 

i 18 1 2 t 4 2293 1 6 3581 1 8 71 1 10 

sino ^ ' ^ ~ ~ 1225 + 8820 ~~ 3922033500 + 1616027212800 ~ 10908183686400 + 
94789 1 12 76477 1 14 

+ 



6178831567324416000 2570452778021883955200 

40722H 16 245265109 1 18 

4 



8412390909889802035200000 365509013631238281 1899727872000 
40956413 1 20 ^ /22 , 

J |_ Q If 22 ) 

510254748374093327017340928000000 v ; ' 

„ N 20 1 2 83 t 4 174931 t G 9605 1 8 

B*. 4;t) = 1 - — + — — - oogoonooi>ygo + 



1323 711480 286339821768 3926946127104 
458505H 10 5892151877 1 12 



572172412994582400 262340410524913476467712 

586063249 t 14 „ , 1fiN 

+ O (t 16 ) , 



10556078423502470838819840 

B . (8!t ) = i-«!f + JZ"L +0 («-). 

S1 „cv ) ; 3267 -r 4008004 v ; 
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